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Abstract 

We provide a geometric method to stabilize asymptotically with phase an ar¬ 
bitrary fixed periodic orbit of a locally generic three-dimensional Hamiltonian dy¬ 
namical system. The main advantage of this method is that one needs not know 
a parameterization of the orbit to be stabilized, but only the values of the Hamil¬ 
tonian and a fixed Casimir (of the Poisson configuration manifold) at that orbit. 

The stabilization procedure is illustrated in the case of the Rikitake model of geo¬ 
magnetic reversal. 
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1 Introduction 

The purpose of this article is to provide a method to stabilize asymptotically with phase 
an arbitrary fixed periodic orbit of a locally generic three-dimensional Hamiltonian dy¬ 
namical system. The main advantage of this method is that one needs not know a 
parameterization of the orbit to be stabilized, but only the values of the Hamiltonian 
and a fixed Casimir (of the Poisson configuration manifold) at that orbit. Moreover, if 
there are many periodic orbits located on the same common level set of the Hamiltonian 
and the Casimir, then the same perturbation can be used in order to asymptotically 
stabilize all of them in the same time. The method can be applied for a large number of 
concrete dynamical systems coming from various sciences, which admit three-dimensional 
Hamiltonian realizations, e.g., Euler’s equations of free rigid body dynamics ([5], [2]), 
the Rikitake system([ll]), the Rossler system ([13]), the Rabinovich system ([14]), etc. 

Before explaining how the stabilization method works, let us clarify what is meant 
by a locally generic three-dimensional Hamiltonian system. Formally, a locally generic 
three-dimensional Hamiltonian system is a three-dimensional Hamiltonian system re¬ 
stricted to an open neighborhood around a regular point of the Poisson configuration 
manifold. Recall that in contrast to Casimir invariants (which globally may not exist), 
locally, around each regular point, every Poisson manifold admits local Casimir invari¬ 
ants (i.e., Casimir invariants of the restricted Poisson structure). Accordingly, a locally 
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generic three-dimensional Hamiltonian system is a dynamical system modeled on a Pois¬ 
son manifold with nontrivial Casimir invariants, each system of this type admitting, 
apart from the Hamiltonian, another constant of motion, namely some hxed Casimir 
invariant (as any two Casimir invariants are functionally dependent, we £x one of them, 
referred to as the Casimir invariant). Due also to the three-dimensionality of the prob¬ 
lem, each orbit of the system is located on a common level set of the Hamiltonian and 
the Casimir invariant. Note that it happens often that there are many obits located on 
the same common level set. If the Hamiltonian or the Casimir invariant is a proper func¬ 
tion, then one expects the dynamical system has periodic orbits. In the case the system 
admits periodic orbits, the stabilization method works as follows: pick a common level 
set containing a periodic orbit of the dynamical system; than the method tells how to 
perturb the system in such a way that the periodic orbit remains a periodic orbit for the 
perturbed system too, the perturbed system keeps dynamically invariant the Hamilto¬ 
nian (or the Casimir invariant) and moreover, the periodic orbit of the perturbed system 
is orbitally phase asymptotically stable with respect to perturbations along the level 
set of the Hamiltonian, which contains the periodic orbit (or along the level set of the 
Casimir invariant, which contains the periodic orbit). Moreover, using a specihc linear 
combination of the above perturbations, one obtains a perturbed system for which the 
periodic orbit is orbitally phase asymptotically stable with respect to perturbations in 
a ’’full” open neighborhood of the orbit. The price to be paid is the loss of the con¬ 
servative nature of the perturbed system. On the other hand, the method can be used 
also in order to obtain the same conclusions, except that the periodic orbit becomes 
unstable for the perturbed system. In all cases, the perturbed families of dynamical 
systems are parameterized by arbitrary strictly positive smooth real functions. This 
method is an improvement of the n-dimensional stabilization technique introduced in 
[9], for the particular case of locally generic three-dimensional Hamiltonian systems. In 
contrast with the general case analyzed in [9], where in order to apply the stabilization 
method, one needs to know a parameterization of the periodic orbit to be stabilized, 
here, the only requirement regarding the periodic orbit is the knowledge of the values of 
the Hamiltonian and the Casimir at the orbit. Moreover, for the class of locally generic 
three-dimensional Hamiltonian systems, we prove that one of the main hypotheses from 
[9] is always fulhlled. 

The structure of the paper is the following. In the second section, we discuss the local 
equivalence between the class of three-dimensional Hamiltonian dynamical systems, and 
that of three-dimensional completely integrable systems. The last section is the main 
part of this paper, here we introduce the stabilization technique and also illustrate the 
main result in the case of the Rikitake model of geomagnetic reversal. 

2 3D Hamiltonian systems versus 3D completely in¬ 
tegrable systems 

In this section we recall the local equivalence between the class of 3D Hamiltonian 
dynamical systems modeled on Poisson manifolds, and that of 3D completely integrable 
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systems. For more details regarding the geometry of Poisson manifolds and related topics 
see, e.g., [3], [5], [1]. 

Let us start by choosing an open set f/ C and a smooth bivector field 11 G = 

T{A'^TU) given locally by 


d d d d d d 

n(x, 2 /,z) = a{x,y,z)—A — + (3{x,y,z)—A — + '-^{x,y,z)—A—, (V)(x,|/,z) e U, 


where a,/3,7 G are smooth real functions. 

Recall that the bivector held 11 is a Poisson structure on U if and only if [11,11] = 0 
(where [•, •] denotes the Schouten-Nijenhuis bracket), this condition being in turn equiv¬ 
alent to the integrability of the smooth 1-form oj = a(ix + (idy + '^dz, namely, dcj Aca = 0. 
In the case when 11 is a Poisson structure on U, the pair {U, 11) is a Poisson manifold. 

Let us recall now the dehnition of the sharp isomorphism between the space of smooth 
1—forms, 11^(17), and the space of smooth vector helds, 3C{U), induced by the canonical 
inner product on More precisely, to each smooth 1—form u = adx -|- (3dy -|- ydz, 
a,/9,7 G we associate a smooth vector held, a;**, given by 


J = a— + (d— + -f—. 
ox oy oz 


In the above notations, given a Poisson structure LI on 17, a three-dimensional dy¬ 
namical system of the type 

h = n(t.)Vi/(t.), 

where H G C°°(17, M) is a smooth real function, and 


n(«) 


0 ~f(u) 

— 7 (m) 0 a(u) 

/3(u) —a{u) 0 


(V)m G 17, 


is called Hamiltonian dynamical system. 

The associated vector held, u E U i— > G TyU = is called the 

Hamiltonian vector field associated to the energy/Hamiltonian H G C°°{U,M.), and is 
denoted by Xh. 

Using the canonical Lie algebra isomorphism between (M^, -1-, -k, x) and (so(3), -|-, -r, [•, •]), 
the so called ”hat” map, ' : —)■ so (3), it follows that 


n(M) = (V)m G U. 


Consequently, using the properties of the ”hat” map we get that for every u E U, 
Xh{u) = f[{u)VH{u) = -(^)VH{u) = -J{u) X VH{u) = VH{u) x J{u). 
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Let us recall now a generic local form of Hamiltonian dynamical systems modeled 
on a three-dimensional Poisson manifold. In order to do that, we first recall a generic 
local form of a general three-dimensional smooth integrable 1—form, uj. More precisely, 
around each regular point (i.e., a point where u in nonzero), u can be written in the form 
UJ := ujy^c = i^dC, for some locally defined smooth real functions C and z/, with v non¬ 
vanishing. The associated Poisson structure will be denoted by As the integrable 

1—form UJ = ydx — xdy do not admits a global representation of the type z/dC, it follows 
that not every three-dimensional Poisson structure can be written globally as Hj.c'. 

Summarizing, each three-dimensional Hamiltonian dynamical system, can be realized 
locally (around regular points) as a Hamiltonian system dehned on a Poisson manifold 
of type (f/, Hj.^c'), where C is an open set, and u,C E M) are smooth real 

functions, with u non-vanishing. More precisely, as (z^dC)** = z/VC, the expression of 
the Hamiltonian vector field generated by a smooth function H G C°°{U, M), is given by 


Xh{u) = S/H{u) X = VH{u) x {udCf{u) = VH{u) x (z/VC)(m) 

= v{u) (ViL(n) X WC{u )), (V)n G f/, 

and consequently the associated Hamiltonian dynamical system becomes 


n?/ 

= v{u) iyHiu) X VC'(n)). 


( 2 . 1 ) 


A straightforward computation shows that H and C are both hrst integrals of the 
Hamiltonian system (2.1). The smooth function C, is called Casimir function of the 
Poisson structure and has the dynamical property that it is a first integral of each 
Hamiltonian system realized on the Poisson manifold (17, Hj^ c”). 

Consequently, any three-dimensional Hamiltonian system, is locally a completely inte¬ 
grable system, i.e., a three-dimensional dynamical system which admits two hrst integrals 
(functionally independent at least on some open subset of the domain of dehnition). Con¬ 
versely, we have that that any three-dimensional completely integrable dynamical system, 
can be expressed at least locally as a Hamiltonian dynamical system of type (2.1). For 
a similar n-dimensional result, see [10]. 

S.ZL 

Indeed, let — = X(m), be a three-dimensional dynamical system generated by a 
at 

smooth vector held X G j£(M^). This system will be completely integrable, on an open 
subset U C if there exist Ii,l 2 G C°°{U,M) two hrst integrals of X, functionally 
independent at least on some open subset of U. 

Hence, the smooth functions Ii,l 2 verify the equalities 


(Cxh) {u) = 0, (Cxh) {u) = 0, (V)n G U, 


or equivalently 


(X(n), V/i(n)) = 0, {X{u), Xhiu)) = 0, (V)n G U. 


( 2 . 2 ) 


From relations (2.2) it follows that for any u E U, the vector X{u) is orthogonal to 
both V/i(m) and Xhi^u), hence X{u) must be collinear with Xlfiu) x Xhiu). 
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Consequently, depending on the existence and the location of the equilibrium points 
of the vector held X, there exists an open set C and a smooth function p G M) 

such that 

X{u)=^Iiu) (V/i(n)x V/ 2 (n)), 

for any u E Q. When dealing with concrete dynamical systems, one often notes that the 
open set is a proper subset of U. 

Consequently, any three-dimensional completely integrable dynamical system is lo¬ 
cally a Hamiltonian system, dehned on a Poisson manifold of the type (f2,nj,c;), with 
C an open subset of M^, z/ = p, C = hln, and energy H = Ii\q. 

Summarizing, in dimension three, completely integrable systems and Hamil¬ 
tonian systems are locally the same. 

Moreover, the generic integral curves of a Hamiltonian dynamical of type (2.1), are 
located on the connected components of the common level surfaces corresponding to 
regular values of H and C. Consequently, if any of the smooth functions C* or if is 
proper, then one expects generic existence of periodic orbits of the system. Recall that 
a smooth function F G is called proper if the preimage of every compact set 

in M is compact in which is equivalent to lim||,^||^oo = oo. 

Note that a large number of concrete dynamical systems from various sciences admits 
three-dimensional Hamiltonian realizations of this type, e.g., Euler’s equations from 
the free rigid body dynamics ([5], [2]), the Lotka-Volterra system([12]), the Rikitake 
system([ll]), the Rossler system ([13]), the Rabinovich system ([14]), etc. 

3 Geometric asymptotic stabilization of periodic or¬ 
bits of three-dimensional Hamiltonian systems 

The purpose of this section is to improve, for locally generic three-dimensional Hamilto¬ 
nian dynamical systems, the stabilization technique introduced in [9], and also to provide 
a new stabilization result. More precisely, we show that one of the main hypothesis re¬ 
quired by the stabilization technique from [9], is always satished in the case of locally 
generic three-dimensional Hamiltonian systems. Moreover, we will show that the elim¬ 
ination of this hypothesis leads to a new stabilization method, which provides phase 
asymptotic stability of a certain periodic orbit, even in the case when we are not able 
to hnd a parameterization of the orbit (which is the case in many concrete dynamical 
systems). The only requirement concerning the periodic orbit to be stabilized is the 
geometric location, more exactly, an implicit characterization of a dynamically invariant 
set containing it. 

Let us start by recalling some dehnitions concerning the stability of the periodic 
orbits of a dynamical system. In order to do that, let x = X{x) be a dynamical system 
generated by a smooth vector held X G X(t/), dehned on an open subset U of an 
n—dimensional smooth Riemannian manifold, whose canonical distance (induced by the 
metric) we denote by dist. Suppose T = { 7 (t) E U :0<f<T}isa T—periodic orbit 
of i = X{x). For more details regarding the stability analysis of periodic orbits see e.g., 
[4], [15], [6]. 
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Definition 3.1 ([9]) 1. The periodic orbit T is called orbitally stable if, given e > 

0 there exists a 6 > 0 such that dist(x(t; Xq), F) < £ for allt > 0 and for all Xq E U 
such that dist(a;o, F) < 5, where 1 1 —)■ x{t; xq) stands for the solution of the Cauchy 
problem x = X{x), a;(0) = xq. 

2. The periodic orbit F is called unstable if it is not orbitally stable. 

3. The periodic orbit F is called orbitally asymptotically stable if it is orbitally 
stable and moreover (by choosing 6 from item (1) smaller if necessary), 

dist(x(t; xo), F) —)■ 0 as t ^ oo. 

4- The periodic orbit F is called orbitally phase asymptotically stable, if it is 
orbitally asymptotically stable and there is a S > 0 such that for each xq E U with 
dist(xo,F) < 6, there exists Oq = do(xo) such that 

lim dist(x(t; xq), 7(t + Of)) = 0. 

t^OO 

We recall now the main theorem from [9] which provides for the class of n - dimen¬ 
sional completely integrable dynamical systems, a constructive perturbation method 
which keeps an a-priori fixed proper subset of first integrals of the system, preserves 
a given periodic orbit, and implements an a-priory prescribed stability behavior of the 
orbit. 

Theorem 3.2 ([9]) Let X = X{x) be a completely integrable dynamical system generated 
by a smooth vector field X G X{U) defined eventually on an open subset U C M”', and 
let k,p E N be two natural numbers, with k + p = n — 1, such that there exist n — 1 
first integrals Ji,..., J^, Di,... ,Dp E C°°{U, M), independent on an open subset V U. 
Suppose the system x = X{x) admits a T—periodic orbit F = {'y{t) E V ■ 0 < t < T} 
such that: 

• F C JD“^({0}), and 0 G is a regular value of the map 

JD = (Ji,..., ...,Dp)-.VCR^^ R-i, 


• VJi(7(t)),..., VJfc( 7 (t)), VF>i(7(t)),..., VF>p(7(t)), X(7(t)) are linearly 
independent for each 0 <t <T. 

If moreover, 0 G zs a regular value of the map .J = (Ji,..., Jk) : Id C M” —)■ M^, 
then for any choice of smooth functions hi,..., hp E M) such that 


/zi( 7 (s))ds < 0, 


hp( 7 (s))ds < 0, 


F, as a periodic orbit of the dynamical system x = X(x) -|- Xo(x), x G V, 


Xn = 


/\ V A A /\ V J, 

i=i 


i=l 


-2 


n—1 




i=l 
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Qi — -k 


p k / p k 

l\ VD,a/\VJzA* /\VD, A/\VJi 

= * = 1 Vj = l l = l 

is orhitally phase asymptotically stable, with respect to perturbations in V, along the 
invariant manifold J^~^({0}). 

On the other hand, for any choice of smooth functions ki, ... ,kp G such 

that there exists io G {1,... ,p} for which 

^io(7(s))ds > 0, 

r, as a periodic orbit of the dynamical system x = X{x) + Xq{x), x eV , 




Xn = 


/\xd,a/\ xjj 


i=l 


i=i 


-2 




n—1 


2=1 


0 ,; — -k 


/\ VA A /\ VJz A ^ VA A A 




1=1 


\j=l 


1=1 


is an unstable periodic orbit. 

Here, * stands for the Hodge star operator for multi-vector fields, and ||Xi A ■ • • A 
Xn-i\\n-i is egual to the Gram determinant of the set of smooth vector fields {Xi ,..., X^-i} C 
X{U). 


Next theorem gives a perturbation of a completely integrable system such that the 
resulting system keeps invariant an a-priori hxed periodic orbit and implements a pre¬ 
scribed stability behavior of the orbit, this time with respect to perturbations located in 
a full open neighborhood of the orbit. The price to be paid is the loss of the conservative 
nature of the perturbed system. More precisely, the perturbed system will not have any 
globally dehned hrst integral, since the characteristic multiplier 1 of the periodic orbit 
to be stabilized has multiplicity one. 


Theorem 3.3 Let x = X{x) be a completely integrable dynamical system generated by 
a smooth vector field X G X(f/) defined eventually on an open subset U ML such that 
there exists a vector type first integral I := (Ji,...,/„_i) G C°°{U,M'^~^) of full rank 
on an open subset V U. Suppose the system x = X{x) admits a T—periodic orbit 
r = { 7 (f) EV :0<t<T} such that: 

• r C and 0 G is a regular value of the map J : 1/ C M"' —)■ 

• V/i( 7 (f)),..., V/n-i( 7 (t)), ^( 7 (t)) are linearly independent for each 0 <t <T. 
Then for any choice of smooth functions hi,, hn-i G C°°{V, M) such that 

[ hi{'y{s))ds < 0 ,..., f A_i( 7 (s))ds < 0, 

Jo Jo 
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r, as a periodic orbit of the dynamical system x = X(x) + Xq{x), x E V, 



n—1 

-2 n-1 

n—1 

/ \ 

Xo = 

Av/, 

■ ^(—where 0* = * 

A A'^^d 


i=l 

n—1 


Vi=i / - 


is orbitally phase asymptotically stable, with respect to perturbations in V. 

On the other hand, for any choice of smooth functions ki,, kn-i E M), such 

that there exists — 1} for which 


^io(7(s))ds > 0, 


Jo 

r, as a periodic orbit of the dynamical system x = X{x) + Xo(a;), x eV , 



n—1 

-2 n-1 

n—1 

/ n—1 \ 

Xo = 

Av/, 

■ where Qi = -k 

A v/,A* /\v/7 


i=l 

n—1 ^—1 

7=1 J A* 

\i=i / . 


is an unstable periodic orbit. 


Proof. As the perturbed vector field Y := X + Xq verifies by construction that 

Cylj = kixlj + ^Xolj = 0 + hjlj = hjlj, (V)j E {1,..., u — 1}, 

the rest of the proof is a direct consequence of the following result from [9]. ■ 

Theorem 3.4 ([9]) Let X = Y{x) be a dynamical system generated by a smooth vector 
field Y E X{U) defined eventually on an open subset U C M”, such that there exist 
k,p E N, k + p = n — 1, and respectively Ji,..., Di,..., Dp, hi, ..., hp E C°°{U, M) 
such that CyJi = ■ ■ ■ = CyJk = 0, and CyDi = hiDi, ..., CyDp = hpDp. Suppose that 
r = { 7 (t) E U : 0 < t < T} is a T—periodic orbit of x = Y{x), such that the following 
conditions hold true: 

• r C and 0 G is a regular value of the map 

JD = {Ji,...,Jk,Di,...,Dp) :U CW^ ^ 


• VJi( 7 (t)),..., VJfc( 7 (t)), VT>i( 7 (t)),..., VT>p( 7 (t)),F( 7 (t)) are linearly 
independent for each 0 <t <T. 


Then, the characteristic multipliers of the periodic orbit T are 


l,...,l,exp 


hi( 7 (s))ds^ ,... ,exp hp{py{s))d^ . 


k-\-l times 














Let us consider now a locally generic three-dimensional Hamiltonian system, i.e., a 
dynamical system of the type 

^ = Kf) X VC(«)), (3.1) 

where G are three given smooth real functions defined on an open 

subset t/ C 

Recall from the previous section that this system is a Hamiltonian system modeled 
on the Poisson manifold (t/, Hj^ c)) and moreover, it is completely integrable, since it 
admits two first integrals, C and H. In order to eliminate trivial cases, we suppose that 
H and C are functionally independent at least on an open subset V O U. 

The following result implies a characterization of the set of equilibrium points of the 
system (3.1), generated by the smooth vector field X{u) = u{u) {'VH{u) x VC(m)), for 
every u E U. 

Proposition 3.5 The vectors'VH{u),'VC{u),X{u), are linearly independent if and only 
if u is not an equilibrium point of the dynamical system (3.1). 

Proof. Note first that the vectors XH{u), XC{u),X{u), are linearly independent if and 
only if 

det{VH{u) I VC(m) | X(m)) ^ 0. 

Using the properties of the triple and respectively the canonical scalar product on 
one obtains that for every u E U the following equalities hold true 

det(Vi/(u) I VC'(u) I X{u)) = det{VH{u) \ VC{u) \ u{u){VH{u) x VC'(m))) 

= z/(m) det{VH{u) I VC'(m) | VH{u) x VC'(m)) 

= iy{u){VH{u), [VC'(u) X {VH{u) x VC'(u))]) 

= iy{u){VH{u) X VC{u),VH{u) x VC'(u)) 

= u{u)\\XH{u) X VC'(m)||2. 

Hence, it follows that the vectors XH{u),XC{u),X{u), are linearly independent if and 
only if i/(m) 7 ^ 0 and XH{u) x VC(m) 7 ^ 0, which is equivalent to X{u) 7 ^ 0, i.e., u is not 
an equilibrium point of (3.1). ■ 

An extremely useful consequence of the above result asserts that the second hypothe¬ 
sis of the Theorem 3.2 is always verihed for locally generic three-dimensional completely 
integrable systems, i.e., for Hamiltonian systems of the type (3.1). 

Remark 3.6 Suppose the dynamical system (3.1) admits a periodic orbit T = { 7 (t) E 
U : 0 < t < T} with period T > 0. Then the vectors XH{'-f{t)),'\/C{'y{t)),X{'y{t)) are 
linearly independent for each 0 < t < T, since periodic orbits contain no equilibrium 
points. 
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Next theorem is the main result of this paper and is an improvement of Theorem 3.2 
and Theorem 3.3, in the case of locally generic three-dimensional Hamiltonian dynamical 
systems. Note that in this case, the second hypothesis of Theorem 3.2 (and Theorem 
3.3) is always verified and moreover, one needs not know an explicit parameterization 
of the periodic orbit to be analyzed. All we need to know is the geometric location of 
the orbit. Since in the case of a locally generic three-dimensional Hamiltonian system 
(3.1), any orbit is located on some common level set of H and respectively C, we shall 
consider two classes of perturbations, namely, one which keeps dynamically invariant 
the Hamiltonian H, and respectively one which keeps dynamically invariant the Casimir 
function C. In both cases, the periodic orbit to be analyzed will remain a periodic orbit 
for the perturbed system too. 


Theorem 3.7 Let 


du 

dt 


= u{u) X VC(m)) , ueU 


(3.2) 


be a locally generic three-dimensional Hamiltonian dynamical system, realized as a Hamil¬ 
tonian system on the Poisson manifold where U C is an open set, and 

H,C,i' E C°°{U, M) are given smooth real functions, such that H and C are functionally 
independent on an open subset V . 

Suppose there exists V (Z V a periodic orbit of (3.2). //T C [H,C)~^{{{h,c)}), 
where {h, c) E M? is a regular value for the map {H, C) : U ^ then the following 
conclusions hold true. 

1. If c is a regular value of the map C : —)■ M, then for every smooth function 

a E C'°°(l/,(0,oo)).- 

(a) T, as a periodic orbit of the perturbed dynamical system 
d?/ 

= u{u) (yH{u) X VC'(m)) - a{u){H{u) - h) [VC'(m) x {VH{u) x VC(m))] , 

Ll6 

u E V, is orbitally phase asymptotically stable, with respect to perturbations 
in V, along the invariant manifold C“^({c}). 

(b) T, as a periodic orbit of the perturbed dynamical system 
(111 

= u{u) (yH{u) X VC'(m)) + a{u){H{u) - h) [VC{u) x {\/H{u) x VC(m))] , 

(-16 

u E V, is unstable. 

2. On the other hand, for every pair of smooth functions a, ft E C°°(l/, (0, oo)).- 

(a) T, as a periodic orbit of the perturbed dynamical system 
du 


dt 


= n{u) {VH{u) X VC'(m)) - a{u){H{u) - h) [VC'(m) x {yH{u) x VC'(m))] 

+ P{u){C{u) - c) [VH{u) X {WHiu) X VC(m))] , 


u E V, is orbitally phase asymptotically stable, with respect to perturbations 
in V. 
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(b) r is unstable, as a periodic orbit of the perturbed dynamical system obtained 
from the above system by changing the sign in front of a or ft. 

Proof. The conclusion follows mainly by Theorem 3.2, Theorem 3.3 and Remark 3.6. 

In order to apply these results, recall hrst from the previous section that the system 

(3.2) is a completely integrable system on the open subset U C since the smooth 
real functions H and C are two hrst integrals of (3.2), functionally independent on an 
open subset V ffU. 

Recall also from Remark 3.6 that WH{u),'S/C{u),X{u) := u{u) x VC(m)) 

are three linearly independent vectors, for every u G T, and hence the second hypothesis 
of Theorem 3.2 is verihed. 

Since H and C are hrst integrals, their values are constants along the solutions of 

(3.2) , and hence as T is a periodic orbit, there exists (h,c) G such that {H,C){u) = 

(h, c), for every u G T. This is equivalent to H{u) — h = 0, and C{u) — c = 0, for every 
M G T, and consequently T will be also a periodic orbit for the above dehned perturbed 
systems, given by items (1) and (2). 

Let us show now that ^“^({c}) is a dynamically invariant set for the perturbed 
dynamical system from (1), i.e., 

QT/ 

= u(u) {VH{u) X VC'(m)) - a{u){H{u) - h) [VC(m) x {VH{u) x VC(m))] G V. 

CJ.6 

This follows if one proves that C is a hrst integral for the above dynamical system. 
Indeed, if we denote for every u E V 

Xo{u) := -a{u){H{u) - h) [VC'(n) x {VHiu) x VC'(m))] , (3.3) 

then the following equalities prove this assertion. 

{Cx+XoC) {u) = {X{u) + Xo(n), VC(n)) 

= (i/(m) {XH{u) X VC(m)) - a{u){H{u) - h) [VC(m) x {XH{u) x VC'(m))] , VC(m)) 

= iy{u){VH{u) X VC(m), VC(m)) - a{u){H{u) - h){VC{u) x {VH{u) x VC{u )), VC(m)) 
= 0 - a{u){H{u) - h){VC{u) x {VH{u) x VC'(m)) , VC'(m)) 

= -a{u){H{u) - h){{VC{u),VC{u))VH{u) - {VC{u),VH{u))VC{u),VC{u)) 

= -a{u){H{u) - h) {\\VC{u)f{VH{u),VC{u)) - {VC{u),VH{u)){VC{u),VC{u))) 

= -a{u)iHiu) - h) {\\XCiu)\mVHiu),XCiu)) - (VC'(n), Vi/(«))||VC'(n)||2) 

= 0. 

Now we have all ingredients necessary to apply the Theorem 3.2 for J := C — c and 
D := H -h. 

(la) In order to prove this item, we apply the Theorem 3.2 for the completely inte¬ 
grable system (3.2), a periodic orbit T C (H,C)~^{{{h,c)}), (where (h, c) is a regular 
value of the map {H,C), and c a regular value for the map C), and the hrst integrals 
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J,D : V —)■ M, J{u) = C{u) — c, D{u) = H{u) — h, where V O U is an open set where 
H and C (and implicitly D and J) are functionally independent. 

If one denotes by X the vector held which generates the system (3.2), then by The¬ 
orem 3.2, we get that for any smooth function ( G C'°°(l/, (0, cxo)), T, as a periodic orbit 
of the system 

3.11 

— = X{u) + Xo{u), uev, 

Xo{u) = \\VD{u) A VJ{u)\\^^ ■ {u))D{u) ^ [VJ{u) A i.{VD{u) A VJ(m))], 

is orbitally phase asymptotically stable, with respect to perturbations in V, along the 
invariant manifold T~^({0}) = (^“^({c}). 

Note that the integral condition from Theorem 3.2 is verihed because Im(—C 
(— 00 , 0 ), and hence one obtains that < 0, for any parameterization 

t G [0,T] t-A 7 (t) G T of the T-periodic orbit T. 

In order to hnish the proof of item (1), it is enough to show that the vector held Xq 
has the same expression as (3.3). 

This follows straightforward taking into account that for any u, u G the following 
equalities hold true: 

u X V = ^{u A u), \\u A v \\2 = II * (u A u)|| = ||m x u||. 

Indeed, since V J = VC and VD = VH, we obtain successively that for any u E V 
Xo{u) = II VD(m) a VJ(m)||^2 • (-1)3-1(-C(m))Z^(m) x [VJ{u) A x{VD{u) A VJ{u))] 

If one denotes 

«('W) := MIO ) U E V, 

^ ^ \\VH{u) X VC'(m)||2’ 

then the expression of vector hied Xq is the same as (3.3). 

Since sign(a) = sign((C), we obtain that {—a{ 7 {t)))dt < 0 for any parameterization 
t G [0, T] I— )■ 7 (t) G T of the T-periodic orbit T, and consequently as all of the hypothesis 
of the hrst part of Theorem 3.2 are verihed, we get the conclusion. 

(16) This item follows by the same arguments as above, except that we replace in 
the expression of the vector held Xq, the control function ( by and then apply the 
second part of Theorem 3.2. 

(2a)/(26) These items follow by applying Theorem 3.3 and using the same arguments 
as those used in order to prove items (la)/(16). ■ 

Remark 3.8 If one realize the system (3.2) as a Hamiltonian system defined on the 
Poisson manifold {U, with Hamiltonian C, then applying the Theorem 3.7 to this 

Hamiltonian realization of the system (3.2), we obtain similar results. More precisely, 
suppose there exists V (ZV a periodic orbit of (3.2). 

//T C {H,C)~^{{{h,c)}), where {h,c) E is a regular value for the map {H,C) : 
6/ —)■ then the following conclusions hold true. 


12 




1. If h is a regular value of the map H : U —)■ M, then for every smooth function 
OL G C°°{y, (0, C)o)); 


(a) r, as a periodic orbit of the perturbed dynamical system 
(i?y 

= „(u) {VH{u) X VC'(m)) + a{u){Ciu) - c) [VH{u) x {VH{u) x VC'(m))] , 

CJ.6 

u E V, is orbitally phase asymptotically stable, with respect to perturbations 
in V, along the invariant manifold H~^{{h}). 

(b) r, as a periodic orbit of the perturbed dynamical system 

(iH 

= u(u) {VH{u) X VC'(m)) - a{u){C{u) - c) [VH{u) x {VH{u) x VC(m))] , 

CJ.6 

u E V, is unstable. 

2. On the other hand, for every pair of smooth functions a, ft E C°°(V, (0, cxo)); 

(a) r, as a periodic orbit of the perturbed dynamical system 
d?/ 

= u{u) (yH{u) X VC'(m)) - a{u){H{u) - h) [VC(m) x {VH{u) x VC(m))] 

CJ.6 

+ l3{u){C{u) - c) [VH{u) X x VC(m))] , 

u E V, is orbitally phase asymptotically stable, with respect to perturbations 
in V. 

(b) r is unstable, as a periodic orbit of the perturbed dynamical system obtained 
from the above system by changing the sign in front of a or ft. 


Next result follows directly from Theorem 3.7 and provides a method to stabilize 
asymptotically with phase an arbitrary hxed periodic orbit of a two-dimensional Hamilto¬ 
nian dynamical system. Before stating the stabilization method, let us recall that a two- 
dimensional Hamiltonian dynamical system generated by a Hamiltonian H E 
dehned on an open set H C is given by 

dn 

— = /i(n)JV'H(n), V E 12, 


where y E and J 


0 1 

-1 0 


Note that any two-dimensional Hamilto¬ 


nian dynamical systems is completely integrable, since the Hamiltonian "H is a hrst 
integral, and conversely, any completely integrable two-dimensional system is locally a 
two-dimensional Hamiltonian dynamical system. 

Let us state now the two-dimensional version of the stabilization method given in 
Theorem 3.7. 
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Proposition 3.9 Let 


(3.4) 


— = fiMS'V'Hiv), V e fl, 
at 

be a two-dimensional Hamiltonian dynamical system, where C zs an open set and 
/i,?{ G are given smooth real functions such that H do not have any critical 

point in some open subset W 

Suppose there exists V <zW a periodic orbit of (3.4). //F C where /z G M 

is a regular value of H, then for every smooth function a G C°°{W, (0, cxo)).- 

1. r, as a periodic orbit of the perturbed dynamical system 

dv 

= y{v)wn{v) - a{v){H{v) - h)VH{v), v eW, 

Ll6 

is orbitally phase asymptotically stable, with respect to perturbations in W. 

2. r, as a periodic orbit of the perturbed dynamical system 

dv 

— = y{v)W'H{v) + a{v){H{v) - h)VH{v), v eW, 

CJ.6 

is unstable. 

Proof. Using the notations of Theorem 3.7, we define ?7 := f2x (—£, e), U ■.= Wx{—e,e) 
(where £ > 0 is a hxed real number), and H{v,z) := T-Liy), C{y,z) := z, n{v,z) := yiy) 
for every {v,z) G U. The rest of the proof follows from Theorem 3.7 applied to the 
periodic orbit T x {0}, and perturbations along W x {0}. ■ 

Let us now illustrate the stabilization result introduced in Theorem 3.7 in the case 
of a Hamiltonian version of the Rikitake system. 

Example 3.10 Let us start by recalling that the Rikitake system is a dynamical system 
which provides a mathematical model for the irregular polarity switching of Earth’s mag¬ 
netic field [8], Since this system exhibits a chaotic behavior, a good knowledge of the 
conservative part, leads to better understanding of the provenience of its complex behav¬ 
ior. Let us consider now the conservative part of the Rikitake system studied in [11], 
namely the dynamical system generated by the vector field 

X{x,y,z) = {yz + (3y) — + {xz - fix) ^ “ ^2/^ e X(M^), 

where fi is a real parameter. 

As proved in [11], the induced dynamical system, 

ii = X{u), u = {x,y, z) & , (3.5) 

admits a three-dimensional Hamiltonian realization of the type (3.2), where 

1 1 

y, z) = 1, H{x, y,z) = - (-x^ + z/^) - fiz, C{x, y,z) = - {x^ + y^) + z^. 
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Note that the maximal set where VH and VC are linearly independent, is the open 
set given by the complement of the set of equilibrium points of (3.5), namely 

V := \ {{(x, 0, /9) : X e M} U {(0, y, —j3) : ?/ G M} U {(0, 0, z) : z G M}}. 

Recall from [11] that there exists an open and dense subset of the image of the map 
{H, C) : ^ R^, such that each fiber of any element {h, c) from this set, corresponds 

to periodic orbits of Euler’s equations. Moreover, any such element is a regular value of 
{H,C), as well as its components for the corresponding maps, H and respectively C. 

Let {h,c) G R^ belongs to the above mention set, and let T C {H,C)~^{{{h,c)}) be a 
periodic orbit of the dynamical system (3.5). 

Then by Theorem 3.7 and the Remark 3.8, the following conclusions hold true. 

1. For every smooth function a G C-(y,(0,cx))).- 

(a) r, as a periodic orbit of the perturbed dynamical system 


x= yz + (3y - a{x,y,z) 
y= xz — I3x — a{x,y, z) 
z= -xy - a{x,y,z) 


- (—x^ + y"^) — (dz — h 


-{-X^+y^)-[lz-h 


(—x^ +y‘^) — jdz — h 


x{-y^ -2z‘^ + 2Pz) 
y (x2 + 2^2 + 2Pz) 

[z(x2 -1/2) _/3(a:2 +1/2)] 


u & V, is orbitally phase asymptotically stable, with respect to perturbations 
in V, along the invariant manifold 

C"^({c}) = G R^ I ^ [x^ + y‘^) = c|. 

(b) r is unstable, as a periodic orbit of the perturbed dynamical system obtained 
from the above system by changing the sign in front of a. 

2. For every smooth function a G C-(y,(0,cx))).- 

(a) r, as a periodic orbit of the perturbed dynamical system 


x= yz + I3y + a{x,y,z) 
i y= xz — fdx + a{x,y, z) 
2 ; = —xy + a(x, y, z) 


(x2 + //2) + 2;2 _ c 
(x2 + 1/2) + ^2 _ ^ 


x\-^ + Dz-D^ 


(x2 + //2) + 2;2 — c 


[—^;(x2 + 1/2) + /5(a;2 — |/2)] 


u G V, is orbitally phase asymptotically stable, with respect to perturbations 
in V, along the invariant manifold 

H~\{h}) = I(x, I/, z) G R^ I ^ (-x^ + y'^) - ^z = h\ . 
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(b) r is unstable, as a periodic orbit of the perturbed dynamical system obtained 
from the above system by changing the sign in front of a. 


3. For every pair of smooth functions a, beC^{V,{0,oo)): 
(a) r, as a periodic orbit of the perturbed dynamical system 

'1 


X = yz + l3y-a{x,y,z) 

'1 


{—x^ +y^) — [dz — h 


X {—y^ — 2z^ + 2(dz) 


+b{x,y,z) 


{x^ + y^) + z^ — c 


y 




y = xz — (dx — a{x, y, z) 

'1 


{—x^ + y‘^) — fdz — h 


y (a;2 + 2z‘^ + 2Pz) 


+b{x,y,z) 


{x^ + y^) + z^ — c 


X 




z = —xy — a[x, y, z 

'1 


(—+ y'^) — (dz — h 


[z{x^ -y^) - /d{x‘^ + y^)] 


+b{x,y,z) 


- + y'^) + z'^ — c - [—z{x'^ + y'^) + (d{x‘^ — 


u G V, is orbitally phase asymptotically stable, with respect to perturbations 
in V. 

(b) r is unstable, as a periodic orbit of the perturbed dynamical system obtained 
from the above system by changing the sign in front of a or b. 
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